Quantum Criticality of ID Attractive Fermi Gas 
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We obtain an analytical equation of state for one-dimensional strongly attractive Fermi gas for 
all parameter regime in current experiments. From the equation of state we derive universal scaling 
functions that control whole thermodynamical properties in quantum critical regimes and illustrate 
physical origin of quantum criticality. It turns out that the critical properties of the system are 
described by these of free fermions and those of mixtures of fermions with mass m and 2m. We 
also show how these critical properties of bulk systems can be revealed from the density profile of 
trapped Fermi gas at finite temperatures and can be used to determine the T — phase boundaries 
without any arbitrariness. 
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I. INTRODUCTION 

Quantum critical phenomena are associated with phase 
transitions at zero temperature as the parameters of the 
system are varied. They are among the most challenging 
problems in condensed matter, since quantum fluctua- 
tions couple strongly with thermal fluctuations in this 
regime. From this viewpoint, ID integrable models ex- 
hibiting phase transitions at T — are particularly valu- 
able, as they can be solved exactly using Bethe Ansatz 
(BA) [l|. These exact solutions will illustrate the mi- 
croscopic origin of their quantum criticality and provide 
unique advance in the study of quantum critical phenom- 
ena and universal scaling theory at quantum criticality. 
On the other hand, recent advances in cold atoms ex- 
periments have provided a highly controlled environment 
to study ID quantum gases in practically all physical 
regimes, thus allowing one to study the predictions of 
Bethe Ansatz. 

Despite much work on ID Fermi gases @, H[, there 
are little studies of quantum criticality using BA solu- 
tions. This is mainly because the thermodynamic prop- 
erties of integrable models at finite temperature are no- 
toriously difficult and present formidable challenges in 
theoretical and mathematical physics. On the other 
hand, the bosonization field theory prediction [4| on the 
low temperature thermodynamics of ID many-body sys- 
tems merely relies on the low-lying excitations which do 
not give proper thermal potentials in quantum critical 
regime. Here, we first present these studies of quantum 
criticality and scaling functions for ID Fermi gas with 
attractive (5-function interaction, which is known to have 
many interesting phases at T = [H-ll|, including the 
so-called Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) [ij} 
like phase. Our purpose is multifold: (I) We present ex- 
act analytic expressions of the pressure, equation of state, 
and other thermodynamic quantities in the strongly in- 
teracting regime which are valid over all the parameter 



regime in current experiments }13j |. These expressions 
will be useful to compare with current results and to ex- 
amine the features of the model confirmed in the experi- 
ment [l3[ . (II) We show that the pressure P takes on an 
intuitive form. It can be regarded as a sum of the pres- 
sures of free fermions and hard core bosons, plus interac- 
tions of clusters of these entities - 2-clusters, 3-clusters, 
etc. (Ill) We show that at criticality, the pressure of 
this system reduces to that of free fermions or mixtures 
of them. The criticality of the FFLO phase when entered 
from different phase is characterized by the latter. (IV) 
We show that using the scaling properties of density of 
each spin component in the quantum critical region, one 
can map out the T = phase diagram of a bulk systems 
using the T > density data of a trapped g pointed 
out recently by one of us [14j. We show how quantum 



criticality can help determine the presence of the FFLO 
like phase, universal Tomonaga-Luttinger liquid (TLL), 
and how the entire T = phase diagram of the bulk 
system can be mapped out from the finite temperature 
non-uniform density profile of different spin components 
in experiments. This new approach to quantum critical- 
ity will open to study of quantum critical phenomena of 
spinor Fermi and Bose gases with higher spin symmetries. 



II. EXACT T = PHASE BOUNDARIES OF 
ATTRACTIVE ID FERMI GAS 

The Hamiltonian of a ID Fermi gas with attractive 
(5-function interaction is K. = H — fiN — hM, 
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N = Eal dx ^a, M = § /dztyfyt " 1^)> where 
a =t,l are the spin labels, c has the dimension of 
(length) -1 . Because of attractive interaction, i.e. c < 0, 
fermions with opposite spin can form a bound pair with 
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FIG. 1: (Color online)Phase diagram in fi — H plane: The 
phase boundaries are represented by solid lines. The V — F, 
V — P boundaries are given in Eq.((5|, and F — PP and P—PP 
boundaries are numerical plot of Eq. Q and ifijl. The blue- 
dashed lines are the extrapolation of the phase boundaries 
Eq. (0 and (JTJl in the strongly coupling regime, which is the 
region below the purple dotted lines, given by y/ jl + 1/2 ~ 1, 

and jl + h/2 ~ 1. Quantum phase transition occurs as the 
driving parameters jl and h cross the phase boundaries at zero 
temperature. The /j, — h phase diagram was found from the 
ground state Fredholm equations in [1, . 



binding energy eb = fi c 2 /4m and spatial extent |c| _1 . 
The parameter 7 = |c|/n, which is the ratio between 
inter-particle spacing (1/n) and width of the bound pair 
|c| , divides the system into a regime of tightly bound 
pairs (7 < —1), and overlapping pairs (— 1 < 7 < 0). For 
systems with spin polarization, both bound pairs and un- 
paired fermions coexist. The equilibration between them, 
even for very small polarization, will make h *~ £;,; and 
we shall later consider temperature regimes, T < £b,h. 
For our discussions, it is convenient to use dimensionless 
quantities where energy and length are measured is in 
units of Eb and c _1 respectively. We shall define jl = n/eb, 
h = H/sb, t = T/eb, n = n/\c\ = 7 -1 , p = P/\ceb\- From 
the BA equation, it can be shown that the strongly cou- 



pling limit 7 >> 1 corresponds to y/H- h/2 << 1 and 

y/jl+ 1/2 « 1. 

ID interacting Fermi gas with contact interaction is 
one of the most important exactly solvable quantum 
many-body systems, and the associating BA solution 
is among the greatest achievements in mathematical 
physics [161 ] . This model contains only one interaction 
parameter c. Yet it is enough to generate a great diversity 
of intricate collective phenomena, including spin-charge 
separation (for c > 0) and oscillatory paring (c < 0). The 
BA solution, however, is formidable. Despite the consid- 
erable analytical effort in deriving the thermodynamic 
Bethe Ansatz (TBA) equations that describe the ther- 
modynamics of the system, they have only been solved 



numerically in most cases 10 L llj, [l3[ and analytically for 
special coupling regimes [z, 15|. 

The T = phase diagram has been worked out by 
Orso [j| and by others 6-8] using BA equations, which 
describe the ground state within a canonical ensemble 
[Ij . Here we study the phase diagram by taking the T — ¥ 
limit of the TBA equations [2J in a grand canonical 
ensemble 17[ and present analytical critical fields in /x — 
H plane. The phase diagram is shown in Figure [T] It 
consists of four phases: vacuum (V), fully paired phase 
(P), ferromagnetic phase (F), and partially paired (PP) 
or (FFLO-\ike) phase. The phase boundaries between 
V — F,V — P, F — PP, and P - PP are denoted as Ha 
to /i C 4 respectively. The quantum phase transitions in the 
ID attractive Fermi gas are determined by the following 
dressed energy equations @, 0] 
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which are obtained from the TBA equations in the limit 
T -> 0. In the above equations a m (X) = ^ ( mc '/2p+\2 • 
The dressed energy e b (A) < (e u (fc) < 0) for |A| < B 
(|fc| < Q) correspond to the occupied states. The positive 
part of e b (e u ) corresponds to the unoccupied states. The 
integration boundaries B and Q characterize the Fermi 
surfaces for bound pairs and unpaired fermions, respec- 
tively. The phase boundary can be worked out by ana- 
lyzing the band fillings with respect to the field H and 
the chemical potential /1 at T = 0. The (V — F) phase 
boundary is determined by the conditions e u (0) < 
and e b (0) > 0, which gives fi c i = —h/2. Whereas the 
(V — P) phase boundary is determined by the conditions 
e"(0) > and e b (0) < 0, which gives /i c2 = -1/2. 

The (F — PP) phase boundary is determined by 
e b (0) < and e u (±Q) = which gives the critical field 
in dimensionless units by 

A*c3 = - 7: - 77- [Q - (2/Uc3 + h + 1) arctan Q] (3) 

2 27T 



with Q = V2^c3 + h [18( , see the solid (F — PP) phase 
line 113 in Fig. Q] The most complicated phase boundary 
indicating the quantum phase transition (P — PP) from 
fully paired phase into the partially polarized phase may 
be determined by the conditions e u (0) < and e h (±B) = 
0, i.e. the Fermi sea of unpaired fermion starts filling. 
Thus we have 
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which provide exact critical field /i C 4, see solid P — PP 
line in Fig. [1] In order to study quantum criticality of 
strongly Fermi gas, close forms of critical fields are essen- 
tial to determine scaling functions of thermodynamical 
properties. They are 
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While Eq.© applies to all regimes, Eq.((6|) and (O are 
expressions in the strongly interacting regime, see the 
dashed lines in Fig. [TJ The above critical fields © and 
([7]) can be also obtained by converting the critical fields 
in the h — n plane, which were found in 0] , into the fi — h 
plane. 



into coupled equations of thermodynamic quantities, ob- 
tained by approximating the dressed energies to the low- 
est order in t, and taking the T\c\ —> oo limit for the 
so-called "string contributions" which describe the effect 
of the spin waves. To describe quantum criticality ac- 
curately, however, higher order terms in t in the dressed 
energy have to be retained. (See footnote below). 
From the TB A equations 0, 0] , it can be shown that the 
dressed energies can be calculated in terms of plolylog 
functions 
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III. EQUATION OF STATE 

The lack of analytic solutions of the TBA equation 
has made calculations of physical properties cumbersome, 
severely limits ones ability to make predictions and to 
identify the physical origin of observed effects. While 
bosonization or Luttingcr liquid theory can provide quali- 
tative information of low temperature properties, they do 
not give equation of states and are accurate only within a 
limited range of temperature. The construction of more 
transparent solutions for thermodynamic functions be- 
comes even more desirable in view of the recent progress 
in the experiments on ID Fermi gas with attractive in- 
teraction, as well as the successes in deducing equation 
of state of Fermi gases from non-uniform density data. 
A close analytic form of thermodynamic functions will 
certainly make comparisons between theory and experi- 
ments easier. 

The thermodynamics of this system has been stud- 
ied analytically p], [U| and numerically 13, 11, HI using 
TBA equations. Ref. Q studies the free energy at very 
low temperatures for the spin balanced case, and shows 
that it is well described by Tomonaga-Luttinger liquid 
theory. The Luttinger description, however, is incapable 
of describing quantum criticality for it does not contain 
the critical fluctuations. TBA equations are a complex 
set of equations in terms of the so-called "dressed ener- 
gies" of bound pairs, unpaired fermions, and spin wave 
bound states. In Ref. 15[, these equations were recasted 



where 
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and Li s (z) — 2~2T=i zk /k s is the polylog function, K — 
i 

-7^/3/2' and J o{ x ) = Efcto (Mp (f ) 2fc comes from the 
so-called "string" or spin wave contributions. The above 
result depends on an important observation that the con- 
volution terms in the TBA equations converge quickly as 
e b (fc) and e u (fc) become greater than zero at low tempera- 
tures. Using the above asymptotic of dressed energies we 
can calculate pressure in a straightforward way. Substi- 
tuting above dressed energies into the pressure per unit 
length p — p b 
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(10) 



and taking integration by part, thus we obtain the fol- 
lowing dimensionless form of the pressure of the system 
as 



p(t,p,,h) =p/(\c\e b ) =p +p 



(11) 



where p b and p u can be interpreted as the pressure of 
the bound pair and unpaired fermions, and are coupled 
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FIG. 2: (Color online) Pressure vs temperature: Comparisons 
between the numerical results from Eq. (|12[) to (|15|l (repre- 
sented by circles) with the analytical result from Eq. (|17l) up 



to the order Y[ u ' b ^ (1st order), Y^"'"' (2nd order) and Y- 
(3rd oder). 

through the following set of equations 
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where the functions F b , f b , and f% are defined as 
f n = Li n (-e^) , ft = Li n (-e^) (16) 



where Vb = 2fl + 1, u u = £l + h/2: Xb, X u are defined as 
in Eq.(]J|) and Eq.^TSJ. 

The derivation of Eq . (fT2")) to (|T5|) is very involved. That 
we present these equations is to prepare for the later 
discussions of the mathematical manipulation needed to 
extract the singularity near the quantum critical point. 
To solve these equations, one substitutes Eqs (|14l) and 
(fT5f into Eq. (fl~2"]) and (p~3|) . This gives two coupled equa- 
tions of p b and p u , which can be solved by iteration. 
From Eqs (|14j) and ()15j) . we can rewrite the functions 
Li3 (— e Xb /*) and Li3 (— e x "'*) in terms of the function 

and f%. After lengthy algebra, the pressures are given 
by p = V2p b + P u with 
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where x — b^u and {^} and {Y^} are given by 
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The expressions of F 3 b and F 3 " are given in Appendix. 

2 2 

The accuracy of Eq. (fTTl) is shown in Figure 2, where the 
result of the expansion is compared with the numerical 
solution of the recast TBA equations. Eq.(fl~2"j) to (fT5|) 
served as equation of state (represented by circles). The 
vertical dotted line represents the Fermi temperature Tp 
in the Rice experiment in ref. 13[, where Tpjtb ~ 0.1. 



Figure 2 shows that for t < 0.1, (corresponds to T < 
Tp in [13]), the pressure is accurately given by the first 
order correction in Eq. dTTl) (to about 1%). At higher 
temperatures, higher order terms are needed. It is worth 
noting that by including terms up to third order, Eq. ([T7|) 
agree essentially with the exact result for t > 0, 5. 
Furthermore, if defining 
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up to the order of Y± , we can rewrite the pressure ()17j) 
as 
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which give insight into understanding cluster-cluster ef- 
fect. From these analytical result of pressures, we see 
that the bound pairs interfere with each other and with 
excess fermions. The F-terms in Eq. (|17[) can be viewed 
as interactions of boson-boson and boson-fermion clus- 
ters of increasingly size. We find that The pressure Eq. 
(JTTj) can reach 95% accuracy for T/e& = 0.2 by including 
upto the terms, and has less than 1% error if one 

includes the Y§ terms, see Figure [2j 

2 

Finally, we note that in the limit |c| — > 00 , the pressure 
(ITU1) gives the following equation of state P(T, /1, H) = 
\ / 2V(2p b ) + V(jiu) with T <C Mi» Mm where 
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is the pressure of a ID Fermi gas with mass m. In this 
limit, for fixed total density n, the effective chemical po- 
tentials for pairs and unpaired fermions are given by 
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Where P = (Nf - N i )/(N t - iVj.) is the polarization. 
This result reflects the fact that the system reduces to 
a mixture of hard core boson (mass 2m) and a gas of 
unpaired fermions (mass m) in this limit. Since hard core 
bosons in ID behaves like fermions, the thermodynamics 
of the system is that of a mixture of two Fermi gases with 
mass 2m and m. 



IV. QUANTUM CRITIC ALITY 
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Quantum critical behaviors are reflected in the sin- 
gularities in thermodynamic quantities such as density 
n = n/\c\, compressibility k = dh/djl, and magnetiza- 
tion M = M/\c\, which are derivatives of pressure p with 
respect to [i and h. Such singularities, however, can- 
not be obtained by directly differentiating the expansion 
Eq. ([T7|) . The reason is that even though Eq.(fT7f gives 
highly accurate value for the pressure with a few terms, 
all higher order terms which contribute insignificantly to 
the pressure have singularities when differentiated with 
respect to fi. To account for these singular structures, 
say, in the density, one must first differentiate Eq. ([T2|) - 
(|15p with respect to fl to obtained coupled equations of 
dp x /dfl and then solves for these quantities by iteration. 
Quantum phase transition occurs as the driving param- 
eters fl and h cross the phase bounaries at zero temper- 
ature. At very low temperatures, T <C £/, (Boltzmann's 
constant ks — 1), the thermodynamics of the ID FFLO 
like phase is governed by the linearly dispersing phonon 
modes, i.e., the long wavelength density fluctuations of 
the two weakly coupled gases. In this low temperature 
regime, spin strings are suppressed. The suppression of 
spin fluctuations leads to a universality class of a two- 
component TLL in the gapless phase of FFLO, where the 
charge bound states form hard-core composite bosons. 
The leading low temperature corrections to the free en- 
ergy give @,[l5 1 

ttT 2 ( 1 1 \ 
/ « /o - — (- + - • ( 24 ) 



Here fo is the ground state energy and Vb,u are the sound 
velocities of pairs and unpaired fermions, i.e. 
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FIG. 3: (Color online) Contour plot entropy S vs chemical 
potential \i from Eq. Eq. (|12|l to (|15|) for different values of 
effective magnetic field [20( . The cross-over temperatures sep- 
arated TLL phases from quantum critical regimes are deter- 
mined from the breakdown of linear temperature-dependent 
entropy. Here TLLf stands for the Tomonaga-Luttinger liq- 
uid of the single component of unpaired fermions. While 
TLLp stands for the TLL of pairs and TLLpp for the two- 
component TLL of FLLO-like states. The left-most dashed 
line separates quasi-classical regime from the quantum critical 
regime. The cross-over temperatures merge at critical points 
as T -> 0. 



Here the Fermi velocity is vp = hnn/m. The TLL main- 
tains below the cross-over temperature at which the re- 
lation of the linear temperature-dependent entropy (or 
specific) heat breaks down. The equation of state Eq. (fT2")) 
to (p~5|) and the TLL nature ([24]) allow ones to quantita- 
tively investigate quantum criticality of the Fermi gas, 
see Fig. [3l 

Using the standard thermodynamic relations, we can 
derive close forms of density, magnetization, compress- 
ibility which allow ones to capture universal low temper- 
ature thermodynamics of the Fermi gas as well as critical 
phenomena. Without losing generality, we can safely ig- 
nore spin wave contribution at quantum criticality due to 
its exponentially small contribution as T — > 0. For our 
convenience in calculating the thermodynamical proper- 
ties, we denote 
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From the equations Eq. (IT2^ - ([T5|) . We derived the total 
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The density in dimensionless scale naturally service 
the dimensionless equation of state which contain 
free-fermion-likc densities with singular behaviour : 
different critical points. The interaction binding en 
rescales temperature. This close form of equation of s 
is very convenient to make fitting of experimental fj 
temperature density profiles of the ID trapped gas within 
local density approximation, see 13, 20j . 
Similarly, magnetization M — (rij- — n\)/2 
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FIG. 4: (Color online) Schematic illustration of universal crit- 
ical phenomena near critical points for the ID strongly at- 
tractive Femi gas, see Eq. (|40[) . The density of bound pairs 
(excess fermions) become a regular part as the background 
meanwhile the one of excess fermions (bound pairs) become 
a singular part. However, there is no back ground density 
for the phase transition from vacuum into the fully-paired or 
fully-polarized phase, i.e. V — > F or V — P. 



bound pairs (excess fermions) become a regular part as 
the background meanwhile the ones of the Fermi liquid 
of excess fermions (bound pairs) become a singular part. 
Near the T — phase boundaries, we find that n and M 
have universal scaling forms 
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and compressibility k = K£b/\c\ with k = dn/d/i 
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can be derived in a systematic way. 

Quantum criticality describes strongly coupled ther- 
mal and quantum fluctuations of matter as quantum 
phase transitions take place at zero temperature. Such 
quantum phase transitions are uniquely characterized by 
the critical exponents depending only on the dimension- 
ality and the symmetry of the system. Here we show 
that scaling functions at quantum criticality can be cal- 
culated from the above close forms of thermodynamical 
properties. Near the critical points, we expand thermo- 
dynamical properties in the limit |/x — fi c \ <C 1. For the 
quantum critical regime, i.e. T > \fi — fj, c \, we find that 
the thermodynamical properties of the Fermi liquid of 



where n 3,n A,a,b,M are constants independent of a 
and t. The expressions are given by 
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In the above equations n D 3 and n 4 are the background 
densities near the critical points [13 and /X4, respectively. 
At quantum criticality, the above densities can be casted 
into a universal scaling form, e.g. l4 l2lH23j . 
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Where the dynamic exponent z = 2 and correlation expo- 
nent 1/ = 1/2 can be read off the scaling functions within 
the expressions Eq.(|3Tj) to (|38|) from which physical ori- 
gin of quantum criticality is conceivable. 

Eq. (j3"Tj) to has a simple interpretation. First of all, 
it is straightforward to work out the critical properties of 
a ID free Fermi gas (single component), as well as for a 
mixture of two different Fermi gases. In the former case, 
the quantum phase transition (as a function of chemical 
potential) is from vacuum to Fermi gas, which we de- 
note as type (a). In the two component case, if the two 
Fermi gases have different particle numbers, then there 
is a transition from vacuum to the majority component 
as a function of Jl (for fixed h), i.e. type (a). There 
is also another transition from the majority component 
to a mixture, which we denote as type (b). Comparing 
Eq. (f3Tj) to (|38|) with the critical properties of Fermi gas 
mixtures, one notes that the transition V—F and V—P is 
that of type (a) where the "fermion mass" for the V — F 
transition is m, and 2m for the V — P transition. In 
the V — P case, this is due to the fact that the tightly 
bound fermion pair acts like a hard core boson, which in 
turn acts as a fermion in ID. In contrast, the transitions 
F — PP and P — PP are of type (b) , where the Fermi gas 
mixture is made of particles with mass m and 2m. This 
critical phenomena is schematically illustrated in Fig . 2) 



t=0.01 
t=0.015 
t=0.02 
- t=0.03 
t=0.04 





-0.45 -0.5 



FIG. 5: (Color online)Magnetization M and density n vs 
chemical potential fi at different temperatures. The inter- 
sections at the right and left panel give the phase boundaries 
of V — P and P — PP transition respectively. 




-0.5 -0.52 -0.54 



FIG. 6: (Color online) Density vs chemical potential: The 
right panel shows the intersection of the density difference 
{n—hw)/\t at different temperatures, which gives the F—PP 



phase boundary. Here, we have uf = — - ,„ Lii (— e 

with A u « v u + £3 fz/z/y/™- ^ becomes the density of a free 
fermion system in the vicinity of F — PP phase boundary, 
i.e. A u reduces to v u . The left panel shows the intersection 
of densities which gives the V — F phase boundary. It also 
shows that the scaling form begins to fail for t > 0.03. 



V. MAPPING OUT THE T = PHASE DIAGRAM 
USING QUANTUM CRITICALITY 



As pointed out in Ref.[14|, the scaling property of the 
density near a quantum critical point enables one to de- 
termine the T = phase boundary of homogenous bulk 
systems from the non-uniform density profile of a trapped 
gas at T > 0. All one needs to do is to plot the den- 
sity profile n(x) of the trapped g function of local 
chemical potential /j,(x) = \x — V(x) at different tem- 
peratures, where V{x) is the trapping potential. At low 
temperatures, these density curves will intersect at the 
same point. The chemical potential (fi c ) associated with 
this intersection point is T — phase boundary. The 
existence of this intersection is due to the fact that the 
singular part of the density is a function of (/i — fi c )/T, as 



seen in Eq. (|3Tj) to Eq. (f38[) . These intersections are seen 
in Fig. [S] to Fig. [SJ which show the numerical solution 
of the recast TBA equations Eq. (fT2"]) to (ITo]) for the n 
and M as a function of [i over a range of chemical poten- 
tial much larger than the scaling region where Eq . (|3 1[) to 
Eq. ([55|) hold. They are equivalent to plotting the exper- 
imental density data in a manner mentioned above. The 
left and right panel of Fig. [5]show the fj dependence of M 
and h across the V — PP phase boundaries (Eq. and 
Eq. ([35|) ) respectively. In Fig. [6j the left panel shows fi 
vs jl across the V — F phase boundary, Eq.(|3Tj). The right 
panel shows h — fiFVsjl across the F — PP phase bound- 
ary, Eq. ([33| . The intersections of these curves yield the 
critical value fi c at T = 0. In these figures, we displayed 
temperature variations from t — 0.01 to t = 0.04. The 
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former corresponds to T/Tp = 0.1 in the recent Rice set 
up, which is the temperature in the current Rice experi- 
ment 
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In general, the slope of densities at the quantum crit- 
icality can be written as the following universal scaling 
form near the critical points: 



(n(t,/i) - n Q )t 



-1/2 



m=Q 



where x 



and the coefficients c m = — f Li i . 



(41) 



The intersection behaviour can accurately determine the 
temperature from the slops near the intersection point. 
For a given polarization (or say fixed polarization), the 
A becomes constant, (see, Eq. ([3TJ) to Eq.(|38j)). Thus the 
slops are uniquely fixed by temperatures. Perhaps this 
scaling behaviour is a good thermometry 24]. 

Finally, we would like to mention that similar plots 
can be constructed for compressibility across all phase 
boundaries, see Fig. [7J For compressibility, it is impor- 
tant (as stressed before) to include all higher order terms 
in Eq. (|17j) . as they are all singular at criticality; and such 
inclusion can be achieved efficiently by taking derivatives 
of Eq (|27[) with respect to jx together with performing 
proper iteration via Eq. ([l"2"T) - (fT5)) . The expression of the 
critical behavior near various phase boundaries are given 
by 




where k Q 3, k Q 4, A4, A5 are given 



^o3 



K 4 



1 x = J_fl+^_^ 

2 f 1 3Vb + 6b\ 



Ab = 



i:\fb 
1 



7T IT I 

2Vb lOfc" 



2V27T 



(46) 



The critical exponents z = 2 and v = 1/2 can be read off 
the universal scaling function 




-0.56 



FIG. 7: (Color online) Compressibility y/tn (dimensionless) 
vs n for h — 1.1 at different values of t. At the critical point 
/Lt C 3, there is a background compressibility. The curves truly 
intersect at a single point after the background is removed. 



VI. CONCLUSION 

We have studied the quantum critical phenomena of 
strongly attractive ID Fermi gas via exact BA solution. 
We have obtained the equation of state with high preci- 
sion from zero temperature up to the temperature scale 
of binding energy. From the equation of state, we have 
obtain the exact scaling form for density, compressibility, 
and spin susceptibility in the vicinity of the T — phase 
boundaries between different phases. These scaling forms 
illustrate the universal TLL signature and physical ori- 
gin of quantum criticality. The excitations near various 
phase boundary is such that their critical behaviors are 
either described by that of free fermions and that mix- 
tures of fermions with mass m and 2m. 

Our exact results can help analyze the recent experi- 
ments on attractive ID spin-imbalanced atomic Fermi gas 
(l3| . It can also help to verify the working of an algorithm 
for determining the T = phase boundary and quantum 
critical behavior based solely on measured density profile 
14j . The BA method proves to be a powerful tool to ob- 
tain exact results for quantum critical behavior. Further 
application of this method to other integrable systems 
including spinor gases may reveal new excitations and 
interesting physics. 
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K ^,T) = Ko + T^-^(^ ¥ ^) (47) 
with an universal scaling function F{x) = Li_i(ic). 
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Appendix: Coefficients 

The pressures ([T21 and (fTlSj) provide the precise equa- 
tion of states for studying two-component Fermi gases 
with population imbalance. By iteration we found the 
analytical equation of state (fT7|) which provides high in- 
sights into many-body effect in the so called FFLO phase. 
The first two coefficients are given in (|18p . The higher 
order correction terms are given by 




which reveal cluster-cluster interacting effect. 



